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Relying on our theoretical approach for the superconducting critical state problem in 3D magnetic 
field configurations, we present an exhaustive analysis of the electrodynamic response for the so- 
called longitudinal transport problem in the slab geometry. A wide set of experimental conditions 
have been considered, including modulation of the applied magnetic field either perpendicular or 
parallel (longitudinal) to the transport current density. The main objective of our work was to 
characterize the role of the macroscopic material law that should properly account for the underlying 
mechanisms of flux cutting and depinning. The intriguing occurrence of negative current patterns 
and the enhancement of the transport current flow along the center of the superconducting sample 
are reproduced as a straightforward consequence of the magnetically induced internal anisotropy. 
Moreover, we show that related to a maximal projection of the current density vector onto the 
local magnetic field, a maximal transport current density occurs somewhere within the sample. The 
elusive measurement of the flux cutting threshold (critical value of such parallel component J c ||) is 
suggested on the basis of local measurements of the transport current density. Finally, we show that 
a high correlation exists between the evolution of the transport current density and the appearance 
of paramagnetic peak structures in terms of the applied longitudinal magnetic field. 

PACS numbers: 74.25.Vs, 74.25. Ha, 41.20.Z, 02.30.Xx 
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I. INTRODUCTION 

The high interest arisen concerning the problem of magnetic flux depinning in type-II superconductors is markedly 
associated with its relevance to technological and industrial applications achieving elevated transport currents with 
no discernible energy dissipation. Thus, the distribution of vortices in a real type-II superconductor is determined 
by the balance between the electromagnetic driving forces perpendicular to the flux lines, and forces pinning the 
vortices to material inhomogeneities (lattice defects, impurities, and finite-size effects}^—. Per unit volume, this reads 
J x B = F p (or J ± B = F p ). 

On the other hand, it is well-known that various striking phenomena occur when a transport current is applied 
to a nonideal type-II superconductor under the presence of a longitudinal magnetic fielder—. First, a remarkable 
enhancement of the critical current density has been observed in a wide number of conventional and high temperature 
superconducting system s 3 ' 7 ' 22 " — . Second, an intriguing negative longitudinal electric field along the direction of the 
transport current has been also reported 8 - - — . Such a resistive structure in the longitudinal field geometry has been 
intuitively understood in terms of helical domains, closely connected to the force-free current parallel to flux-lines 
(recall that the magnetic force per unit volume is given by J x B) 11 ' 18 ' 21 ' 25 . 

Nevertheless, although several facts suggest that the vortex lattice is arranged in a helical configuration, perhaps 
close to the force-free arrangement (J||B), the arising voltage cannot be straightforwardly described as a critical 
flux-flow voltage if one is to be consistent with the experimental report s 13 ' 32 . On the other hand, although several 
proposals have been done in terms of the crossing and recombination of adjacent nonparalcll flux lines, the so-called 
flux line cutting phenomenon has been basically recognized as the physical mechanism by which the longitudinal 
voltage is produced ^ 20 ' 30 ' 31 ' 33 " — . Outstandingly, remarkable numerical and conceptual difficulties related to the 
implementation of the above picture of the local dynamics for the transport current in realistic longitudinal geometries 
have lead to unfortunate omissions of the phenomenon in practical calculations related to the applications of type-II 
materials. 

Recently, we have shown that our general critical state theory (Refi^), is able to describe on a quantitative 
basis the anomalous features involved in the local electrodynamics of a current carrying superconducting sample 
subjected to variations of both longitudinal and transverse magnetic fields (Ref.— ). In brief, our work is developed at a 
phcnomcnological level which allows to deal with 3D-magnetic field configurations and nontrivial constraint conditions 
defining anisotropy effects on the superconducting material law. We have introduced a geometrical description for 
the constraints on the critical current density in terms of a closed region A r , such that the physically admissible 
states are given by the condition J £ A r . When the physical limitations concerning the microscopic phenomena of 
flux depinning and cutting are considered to be independent, the region A r becomes a cylinder, with a longitudinal 
rectangular section of size 2 J c ± x 2 J c || . This ansat££- can be identified with the so-called double critical state model in 
type-II superconductivity. Hereafter, a related anisotropy parameter, that we will name after bandwidth x = J c \\/JcX 
will be used. 

Preliminar work on the intriguing effects associated to the longitudinal transport problem was presented in Ref£. 
Here, supported by numerical simulations that cover an extensive set of experimental conditions, we put forward a 
much more complete physical scenario. Thus, we will show that the striking existence of negative flow domains, local 
and global paramagnetic structures, emergence of peak-like structures in both the critical current density and the 
longitudinal magnetic moment, as well as the compression of the transport current in type-II superconductors under 
parallel magnetic fields, are all predicted by our general critical state model. In addition, we shall introduce some 
ideas that could be applied for the determination of the flux cutting threshold from local measurements of the current 
density flowing along specific layers of the superconducting sample, as correlated to the behavior of the magnetic 
moment components. 

The paper is organized as follows. In Sec. UH the physical background for the critical state concept in longitudinal 
geometries under applied 3D magnetic fields is introduced (Sec. Ill A|) . Then, the underlying approximations and the 
variational statement for the longitudinal transport problem in superconducting slabs are both described in detail 
(Sec. Ill Bp . This encloses the time-discretized description of the magnetic field penetration process, and the local 
dynamics of the transport current in an arbitrary magnetization process. The theory is then applied for a set of 
experimental configurations and a collection of material laws as regards the interplay of the flux deppining and 
cutting mechanisms (Section IIII[) . On the one hand, the infinite bandwidth model (or T-state model) with J c m 3> J c ± 
is assumed (Sec. IIII A[) . As the necessity of including a model with a well defined value for the threshold J c m will 
become apparent, in sections IIII Bl and IIII CI a more general description is presented as function of the anisotropy 
effects on the material law. These effects will be investigated in terms of the range of parameter values: % = 1, 2, 3, 4. 
The procedure will reveal both local and global properties of the magnetic moment and the transport current flow 
indicating a possible reconstruction scheme of the underlying material law. Sec. HVl is devoted to summarize the main 
findings of our work. 
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II. LONGITUDINAL TRANSPORT PROBLEM IN THE DOUBLE CRITICAL STATE REGIME 

This section is devoted to introduce the theoretical background that justifies our variational statement as an 
appropriate tool which allows to implement a wide class of material laws in the investigation of the electromagnetic 
response of superconducting samples subjected to 3D magnetic field configurations and transport currents. The 
characteristic mathematical equations for the critical state problem are introduced in terms of a variational statement 
that is applied within the infinite slab geometry. The relevance of the parallel and perpendicular critical current 
limitations will be demonstrated. 



A. General statements 



The dynamics of the local magnetic field and transport current profiles in type-II superconductors is commonly 
obtained within the critical state model framework. In the basic formulation (Bean's modelj^^r— ) one considers 
that the magnetic flux lines, when driven by a macroscopic current density, will either penetrate or exit across the 
superconducting surface. In more detail, when the driving forces (local value of J x B) overcome the pinning force 
in a local section of the superconducting sample, the system of vortices rearranges itself into a new metastable state 
such that the vortex lattice is pinned again when the equilibrium at the boundary is reestablished. Notice that this 
physical mechanism relates to the component of J perpendicular to B. Since the displacement of the flux-lines takes 
place with high associated resistivity, the system quickly adjusts itself to successive equilibrium states so as to avoid 
the high resistive losses. On the other hand, reversible energy terms related to the equilibrium properties of the vortex 
lattice are usually neglected (justified by the range of interest for the local magnetic fields H c \ <C H <C -ff C 2), and the 
linear relation B = y«oH is assumed. 

For the particular case that the transport current is applied along the direction of the vortex lattice, rotation of 
vortices is induced by the component of the current density parallel to the magnetic field. Subsequent cross-joining 
and cutting will occur if the angle between neighboring vortices overcomes some critical value. Thus, in general, for 
nontrivial configurations, and related to the threshold values for the current density that trigger the displacement of 
vortex lines in some sense, a critical state may be defined as a given configuration which is able to withstand a critical 
current defined by J c = J c m + 3 C ±, as long as neither threshold (J c m or J c j_) is exceeded 4 -"— i 30 i 31 i 35 i 39 . 

One must mention that, data from a wide number of experiments have shown that in many cases J c m and J c j_ 
can be of the same order of magnitude£ ~ 10 ' 12 i 13 ' 24 i 40 ~ — . Then, in order to explain the accompanying phenomena, 
both thresholds have to be included in the theory. In these cases, the local current voltage characteristic, connecting 
the local electric field E and current density J, is strongly anisotropic with respect to the direction of the magnetic 
induction B. This fact introduces strong difficulties when trying to find an analytical solution for the critical- state 
equation V x H = J c (or 0). Our formulation can somehow bypass such difficulties by treating the electric field as a 
derived quantity that may be obtained a posteriori. 

Thus, the critical-state equation can be posed in numerical terms, assuming an evolutionary discretization scheme. 
Let us suppose that H; stands for the local magnetic field intensity at the time layer I6t, and that the current density 
profiles relate to some magnetic diffusion process that takes place when the local conditions for the material law 
J\\ 5= «7 c ii , Jj_ < J c ± are violated. The successive field penetration profiles within the superconductor may be obtained 
by the finite-difference expression of Faraday's law, 

Mo Hi+1 ~ H ' = -V x E = -V x [p(V x H, +1 )] . (1) 
or 

Here, p(3) plays the role of a nonlinear and nonscalar resistivity that should properly incorporate the physics of the 
threshold and dissipation mechanisms mentioned above. Notice that the local profile H; + i can be solved in terms of 
the previous field distribution Hi and the boundary conditions at time layer (I + l)6t. The initial condition fulfills 
Ampere's law V x H; = J/ as well as V ■ H; = and V • J; = 0. One possibility for making the integration of this 
system affordable is to find an equivalent variational statement. Thus, one can avoid the integration of these set of 
differential equations by inversion of the Euler-Lagrange equations 

J ;+1 -VxH ;+1 =0, (2) 

and 

Vxp ; +H ;+1 -H,=0, (3) 

for arbitrary variations of the time-discretized local magnetic field, H/, that plays the role of the unknown. On the 
other hand, we have introduced the Lagrange multiplier, p ; , which can basically identified with the electric field of 
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the problem. The related Lagrange density (over whole space) is 

L = i|H, +1 - H ; | 2 + p • (V x H /+1 - J m ) . (4) 

As indicated above, minimization under the constraint V x H; + i = J;+i (Ampere's law) allows to reduce the number 
of variables in the statement of the problem, skipping the explicit inclusion of E. 

For practical purposes, we stress that the above equations may be transformed so as to get a vector potential 
formulation by means the condition 

H i+1 H ; = Mo X V x (A m - A,) . (5) 

Thus, by using standard electromagnetic manipulations, one may show that the functional to be minimized can take 
the following form in terms of the current density 

T [J I+1 e A r ] =— / AA • 3 l+1 dr 
Mo Jq 

Ji+1 ', [J,+1 ~ 2J,] dr^ 
oxn l r — r I 

(6) 

It must be emphasized that Eq.© can be applied for any shape of the superconducting volume fl as well as for 
any general restriction (material law) for the current density Jz+i G A r . It is also to be noticed that solcnoidality 
(V • J = 0) has to be imposed so as to be consistent with charge conservation in quasi-steady regime. 



B. Variational statement in slab geometry 



In this subsection, we derive a specific variational formulation for the longitudinal transport problem in super- 
conducting slabs for eventual 3D-field configurations (Fig. [T|), i.e., both in plane and transverse local magnetic field 
components emerge as derived effects of the external sources (transport currents and/or external magnetic sources). 

We shall consider the time evolution of magnetic profiles H/ + i(z) within an infinite superconducting slab (of 
thickness 2a), cooled under the assumption of an initial state defined by a uniform vortex lattice perpendicular to 
the external surfaces, i.e., a constant magnetic field H z q. Then, the constraint for the current density (material 
law J S A r ) corresponding to the limitations on the parallel and perpendicular directions mentioned above may be 
visualized by a finite cylinder with its axis parallel to the local magnetic field H z (Fig.QJi). When a transport current 
is injected along the superconducting slab in the direction of y-axis, a magnetic field component H x appears, inducing 
a rotation of the critical current region A r (see Fig. [!})). Finally, if a third magnetic field component H y is switched, 
a new rotation of the region occurs (Fig. QJ). Notice that, by symmetry, the current density is confined to the XY 
plane, i.e.: J = ( J x (z), J v (z), 0). In particular, this means that in practice one should impose the restriction that J 
belongs to the projection of the critical current region onto the plane (J £ A p ) and that V • J = 0. 

At this point, we must mention that for numerical convenience, the material law is not strictly used in the form of 
a cylinder, but smoothly reshaped as a superellipsoid^^ by means of the relation 

(4r + (£T si - 

The critical current region A r will be characterized in terms of its bandwidth \ = |Jc||/«^c±| an d the superelliptic 
index n. 

We call the readers' attention that two analytical approaches for the slab geometry may be found in the literature 
for extreme situations. The first one was introduced by Brandt and Mikitik in Ref4i for the regime of strong pinning 
with very weak longitudinal current conditions, i.e., H z must be very high as compared to H xy (a) (then Jm <C J c \\)- 
On the other hand, the opposite limit (H z —> 0) was recently considered in Ref.— . Here, and based on the numerical 
resolution of the variational statement, a complete tour along the whole set of values for the perpendicular field will 
be presented. 

In order to simplify the mathematical statements we shall normalize the electrodynamic quantities by defining 
h ee H/J c j_a, j = J/J c j_, and z ee z/a. In turn, our problem will be described in terms of N s discretized layers, each 
one characterized by a current density function j(zj) = J x (zi) + j y (zi) distributed along \zi\ < Nia/N s . For further 
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consideration, notice that the value of the transport current assumed in this work, i.e.: penetration to half thickness, 
imposes the boundary condition h x (a) = 0.5 (Fig. [2]). 

In our discretized description, and in terms of Ampere's law h x (z{) will be evaluated from 



2j»(zj) -jy(zi) 
j<i 

Similarly, the local profiles for the longitudinal magnetic held component h y (zi) can be obtained from 

2 Jx(Zj) ~ jx{zi) 

j>i 



h y (zi) = 8^ , (9) 



with 6 = a/N s the thickness of each layer over the local plane Zi- 

Then, for the magnetic process sketched in Fig. [TJ;, Eq.© takes the following form in terms of discrete variables 



m = i + 1)= \ E I ?,i+iMr j ii l+1 E ZiMijfy+i 

+E 7 u + i(* - V2)(ftr+i(«) - W) ■ ( io ) 



In this expression we have introduced the sheet currents Ifi + i = dj x (zi,t = 1 + 1) and /f ;+1 = 5j y (zi,t = 1 + 1), 
and their mutual inductance coefficients M?- . J- has to be minimized along the different time steps (I = 1,2,...) 
in the magnetic process under consideration. On the other hand, one can show that the parallel and perpendicular 
projections of the sheet current components are given by 



I\ = (1 - //;,)/;., + (l - Irjll - /,.,/.,. 

"x,i ± x,i "T" "y,i ± y,i ~r ^hx,ihy,i ±x,i±y,i 



I\\ — h x i I xi + h i Iy i + 2h Xt ih y> i I x ,il v ,i ■ (11) 



These components have to be constrained according to the material law in Eq.([7]). Further, for the transport problem, 
one has to consider the external constraint: 

XX* = J *-- ( 12 ) 

i 

The theoretical framework becomes closed by the following expressions for the mutual inductance coupling 
elements^ 

Mfj =l + 2[min{i,j}} , 

Mfj = 1 + 2 [N s - max{i, j}] , 

JW* = 2 (i+JV a -i) . (13) 

In summary, the magnetic response of the superconductor is characterized by a collection of discretized current 
elements for the planar sheets (j X i,j V i) at the time steps I + 1 = 1, 2, 3, .... Then, the magnetic field profiles may be 
obtained by means of Eqs.© and ©, and the sample's magnetic moment per unit area from 

M /+1 =Ezxj,, (14) 



where a factor of 2 has been introduced related to the contribution of the U-turns at infinity for each circuit. 
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FIG. 1. (Color online) Scheme of the time dependence of the magnetic process considered in this work, (a) Left panel: an 
external magnetic field H z o is applied perpendicular to the external surface of a superconducting slab of thickness 2a. (b) 
Middle: the slab is later subjected to a transport current along the y-axis. This generates a local magnetic field component 
H x . (c) Right panel: at t = t" an external magnetic field H y parallel to the transport of current direction is switched on. 
The lower row of each panel depicts the region A that constrains the flow of electrical current, as related to the critical state 
material law. A is a cylinder with axis parallel to the local magnetic field, with radius J c ± and length 2J u. The components 
of J within the sample's plane are indicated for each case 



III. INFLUENCE OF THE MAGNETIC ANISOTROPY IN THE CRITICAL CURRENT 



In the previous sections we have visualized the double critical state model as a phenomenological approach which 
can be formulated by means of a variational problem with physical constraints. Here, based on the above mentioned 
theoretical statements for the longitudinal transport current problem, we will show the theoretical predictions for 
the magnetization process outlined in Fig. QJ as h y (a) = H y o is increased. In addition, several initial states h Z Q will 
be focused on. Moreover, we shall concentrate on the effect of the flux cutting boundary (j c \\) considering several 
conditions for the material law. Two extreme cases (x = 1 and \ °°) will be considered first (Sections IIII Al and 
IIIIB[) . and the range in between at a second step (x = 2, 3, and 4 in Section UlI C[) . Remarkably, our procedure will 
reveal the fingerprints of the cutting and depinning mechanism, thus being a theoretical pathway for the reconstruction 
of the material law, represented by the proper region A r . 

Henceforth, we shall use the simplified notation T or CTx as regards to the infinite bandwidth model (T by 
transport) or double critical state model (CT by cutting and transport) with anisotropy x = IJc||/Jc_l|- 



A. T-states for the longitudinal transport configurations 
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FIG. 2. (Color online) Profiles of the magnetic field components h x [z,h y (a)] and h y [z,h y (a)], and the corresponding current- 
density profiles j y [z, h y (a)] and j x [z, h y (a)] for the T-state model and perpendicular magnetic field components h z o = 10 (top), 
hzo — 2 (middle) and h z o = 0.5 (bottom). The different curves correspond to the following sets of values for the longitudinal 
field at the surface: (i) top row: h y (a) = 0.005, 0.050, 0.170, 0.340, 0.500, 0.680, 0.845, 1.0, 40.0, 80.0, 150.0, 300.0, 500.0, 750.0, 
1000, (ii) middle row: h y (a) = 0.0050.050,0.170,0.340,0.500,0.680,0.845,1.0,5.0,10.0, 20.0,30.0,40.0,50.0,60.0, (iii) bottom 
row: h v (a) = 0.001, 0.050, 0.170, 0.340, 0.500, 0.680, 0.845, 1.0, 2.0, 3.0, 5.0, 7.0, 8.0, 9.0, 10.0. 



1. Field and current density penetration profiles 

Fig.|2]shows the magnetic field profiles and the induced currents for three different initial conditions, i.e., h z o = 10, 
h z o = 2, and h z o = 0.5, all of them under assumption of the T-state model. The initial state for the transport current 
condition (I tr = J c ±a/2) establishes the initial transport profile j y {0 < z < a/2} = and j y {a/2 < z < a] = 1. As 
the transport current is no longer modified, the condition h x (a) = 0.5 can be applied in what follows. On the other 
hand, by symmetry, one has to impose the condition h x (0) = at the center of the slab. 

Then, when the external magnetic field h y (a) is linearly increased from h y (a) = 0, a current density j x is induced 
from the superconducting surface as an effect of the Faraday's law. Simultaneously, the local component of the 
magnetic field h x (z) increases monotonically following two continuous stages fulfilling the aforementioned boundary 
conditions. First, the superconducting sample is fully penetrated by the transport current when h*(a) = 0.845 ±0.003 
and eventually, the condition j y (0) = 1 is reached as soon as h y {a) — > 0.860. We notice that the value of h*(a) 
for the full penetration profile is basically independent of h z o, at least to the numerical precision of our numerical 
calculations. This agrees with the analytical solution of Ref.— . Second, a remarkable enhancement of the transport 
current density occurs around the center of slab as h y (a) increases over h*(a). Furthermore, an eventual negative 
current density appears shielding the positive transport current around the center of slab. In more detail, notice that 
the appearance of negative current flow is enhanced when the magnetic component h Z Q is decreased (Fig. [3]). 

Outstandingly, profiles of magnetic field reentry (paramagnetism in the component h y around center of slab) are 
obtained for h Z Q < 1 under relatively low applied magnetic fields h y (a) (see Figj2|). 

Another remarkable property is that, for the range of values h z (0) < h* y (a) negative surface current appears even 
for the partial penetration regime, e.g., for h z o — 0.5 one has j y (a) < for h y (a) > 0.722. Recall that, in Ref£ we 
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FIG. 3. (Color online) Evolution of the local current density as a function of the applied longitudinal magnetic field h y o — h y (a) 
along the central and external sheets of the slab. The results are shown for the T-state model {J c \\ — > oo & J c ± = 1.0). Top: 
the components j y and j a at (z = 0) and (z = a). Middle: details of the above behavior. Bottom: behavior of the parallel 
and perpendicular components of j in the same conditions as above. The different curves correspond to the values of the 
perpendicular magnetic field given by h-o = 1, 2, 5, 10, 20, 50, 100, 200 (all plots having the same color scale). 



have analytically shown that both effects, local paramagnetism and negative current zones are also predicted in the 
limiting case h z o = . Along this line, as a general rule, we can conclude that the smaller the value of h z o, the sooner 
the surface of negative transport current and even paramagnetic local effects appear. 



2. Analysis of the current density behavior 

Fig.[3]displays the evolution of the current density vector as a function of the longitudinal magnetic field h y (a). Note 
that we focus on the specific values at the superconducting surface (z — a) and at the center of the superconducting 
slab (z = 0). 

In order to understand the physical mechanisms responsible for the observed behavior, it will be useful to consider 
the following representation of the current density components. We identify the following decomposition of the vector 
J 

J = J || +J ±(1 +J_l , (15) 

that assumes a polar axes representation, with the parallel, azimuth and polar components of J defined in terms of 
the magnetic field direction. The following expressions are obtained for such components, when a cartesian coordinate 
system is used. 

(i) The current density parallel to H or so-called cutting current component Ji | : 

J]] = H x J x +H y J y (ig) 
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FIG. 4. (Color online) The magnetic moment (M x ,M y ) of the slab as a function of the applied magnetic field com- 
ponent h y (a) for the T-state model. The curves are labelled according to the perpendicular magnetic field component 
h z o = 1, 2, 5, 10, 20, 50, 100, 200, the same color scale applies to all the plots. 



(ii) The component of J perpendicular to the plane defined by the vectors z and H or so-called azimuthal current 
component 

J —HyJ X + H X Jy 

J±a - (fl2 + flj)V> ' [U) 

(iii) the component of J perpendicular to H and contained in the plane defined by the vectors z and H or so-called 
polar current component J±$: 

J _ H Z (H X J X + Hyjy) 

J±e H{H* + H$)W ■ (18) 

Within the T-state model, the behavior of J is straightforwardly explained: the unbounded parallel current density 
allows unconstrained rotations for the flux lines as the applied magnetic field increases. In particular, this leads to 
negative values of j y (a) (slope of h x (a)), simultaneous to high j y (0) (slope of h x (0)). We call the readers' attention that 
negative values of the transport current are favored by smaller and smaller values of the field component perpendicular 
to the surface of the sample h z o- 

A property to be noticed in Fig. |3] is that , at the center of slab the flux line dynamics is mainly governed by the 
longitudinal transport current density j y (0). The basic idea is that for moderate values of h z , when h y increases j y 
practically becomes jn . As this component is unconstrained, it grows indefinitely at the center. 



3. Magnetic moment of the sample 



Let us now concentrate on the magnetostatic properties by means of the global sample's magnetization curve M(H). 
Thus, we have calculated M as a function of the longitudinal magnetic field h y (a). Fig. |4] displays the magnetic 
moment components M x (h y o) and M y (h y o) in units of J c ±a 2 . Notice first that within the partial penetration regime 
(h y (a) ^ h*(a)) the magnetic moment components are almost independent of the transversal magnetic field h z o (at 
least for non small values of this quantity). On the contrary, when h Z Q < 1 and the patterns of negative current 
even occur before of the full penetrated state, magnetization slightly increases. This is accompanied by faint field 
reentry effects that are also shown in the figure. Furthermore, as the threshold cutting current j c u is unbounded for 
the T-state model, the magnetic moment M x always increases as related to the diverging behavior of j y (0). 
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We must emphasize that the unbounded behavior for the parallel current density assumed above (that leads to the 
prediction of arbitrarily high values of the transport current density) must be physically reconsidered. Thus, the trend 
of the magnetic moment M x and also the unbounded longitudinal current density j y disagree with the experimental 
evidences recollected in Refs£ ~ 10 i 18 ' 21 i 25 ~ — . In addition, in Fig. [3] one can notice that, as soon as the flow of negative 
current along the superconducting surface is reached, it never disappears notwithstanding the longitudinal magnetic 
field remains increasing. By contrast, the disappearance of the patterns of superficial negative current were detected 
in Refs i 11 ! 14 ' 15 ! 18 ' 21 . These observations have lead to consider J c ^ -bounded descriptions as satisfactory solutions of the 
peculiar phenomena involved on the longitudinal transport current problem 8 - ^ 11 ' 15 . Rather recent experimental data 
on high temperature superconductors^ also indicate that physical bounds are to be considered for both components 
of the critical current. 

On the other hand, the method described in this work suits the necessity of dealing with a physically acceptable 
description of both local and global issues about the electromagnetic quantities involved on the longitudinal transport 
current problem^. More realistic models for the material law are presented below. 



B. CT1 states as a general approach 



In this section, we show the results obtained for the square condition given by \ = 3c\\/icX = 1 (CT1 in what 
follows). This can be considered as a lower bound for such quantity because the experimental values reported in the 
literature arc typically above unity. 

In order to obtain continuity with the T-state results obtained above, the electrodynamic quantities of interest have 
been proceeded under the same arguments developed in Sec. Ill Bi as regards to the magnetic process shown at Fig. [TJ;. 
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FIG. 5. (Color online) Profiles of the magnetic field components h x [z, h y {a)\ and h y [z, h y (a)] for a perpendicular field component 
hzo — 0.5. Also included are the corresponding current-density profiles j y [z, h y (a)] and jx[z, h y {a)] for the CT1 model. For clar- 
ity, the longitudinal magnetic field is applied in three stages: (Top) h y (a) = 0.005, 0.050, 0.170, 0.340, 0.500, 0.680, 0.845, 1.0, 1.1, 
1.3, (Middle) h y (a) = 1.3, 1.6, 1.9, 2.2, 3.0, 4.0, 5.0, 6.0, and (Bottom) h v {a) = 10, 20, 40, 80, 100, 125, 150, 1000. 
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FIG. 6. (Color online) Same as Fig[5] but for h z o = 10 and the values of the longitudinal field: (Top) h y (a) — 
0.005, 0.050, 0.170, 0.340, 0.500, 0.680, 0.845, 1.0, 4.0, 6.0, 8.0, 11, (Middle) h y (a) = 13, 15, 17, 20, 25, 35, 50, 65, (Bottom) h y (a) = 
80, 100, 125, 150, 200, 300, 400, 1000. 



On the other hand, with the aim of getting a detailed physical interpretation on how the longitudinal and transverse 
magnetic fields affect the dynamics of the transport current problem, we show the magnetic penetration profiles for 
low and high perpendicular fields, i.e., h z o = 0.5 (Fig. [5]) and h z o = 10 (Fig. [6]). Eventually, for completeness, the set 
of initial conditions h z o is extended in (Fig. [7]) . It will be shown that the fingerprint of the CT model is identified as a 
peak effect in the magnetization curves (Fig. [5]) caused by the maximal enhancement of the critical current transport 
density along the central layer. 

In order to ease the interpretation of the intricate behavior of the magnetic profiles under CT conditions, we have 
divided the experimental process in three successive stages as the longitudinal magnetic field component h y (a) is 
increased, i.e.: 

(i) The current density at the center j y (0) increases until a maximum value is obtained (top of Figs. [5] and [5]). 
Notice that the partial penetration regime is included within this stage . 

(ii) The minimum value for the longitudinal current density along the superconducting surface is reached (middle 
row of Figs. [5] and 

(Hi) The longitudinal current density j y (a) stabilizes around j y (a) ~ 0.5 (Bottom of Figs. [5] and [5]). 



1. Field and current density penetration profiles 



To start with, notice that the trend of the profiles for the partial penetration regime is quite independent of the 
perpendicular magnetic field h z o (Figs. [5] and [6]). Moreover, the partial penetration regime in which the transport 
current zone progressively penetrates the sample is practically independent of the magnetic anisotropy of the critical 
state (compare to Fig. [2]). On the other hand, the negative current patterns are also found under the low applied 
magnetic fields h z o < 0.5. However, by contrast to the results within the previous section, recall that for the T-state 
model the condition j c u — > oo allows unbounded values for the longitudinal current j y at the center of the sample. By 
contrast, for the bounded case, the magnetic anisotropy of the material law A r defines the maximal current density 
for the critical state regime. In other words, the maximal length of the vector J within the region A r defines the 
maximal transport current allowed in the superconductor. 
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FIG. 7. (Color online) Dynamics of the current density vector for the CT1 model. To the left, we show the patterns of transport 
current over the central layer (j y (z = 0)) and external layer (j y (z — a)) of the slab. To the right, we show the variation of 
the perpendicular component j x (z = 0) and j x (z = a). The curves are labelled according to the perpendicular magnetic field 
component h z0 = 0.1, 0.5, 1, 2, 5, 10, 20, 50, 100, 200. 



Recalling Refj^£, the maximal longitudinal transport current density corresponds to the optimal orientation in 
which the biggest distance within the superelliptic region is obtained, i.e.: 

/ ,, n\ (n— l)/2n 

jrnax = (l + x 2„/(„-l) j _ (lg) 

Note that his equation allows to obtain the maximum value expected for j y in terms of the actual critical state model 
in use. Thus, for \ = 1, n = 4 one gets j™ ax = 1.3, a value that can be checked in Figs. [5] and[6] 



2. Analysis of the current density behavior 



Complimentary to the field penetration profiles, Fig. [7J shows the dynamics of the main components of the current 
density along the central and external layers of the superconducting slab in the CT1 condition. Once again, notice 
that the full penetration regime can be clearly distinguished from the partial penetration regime. 

Also, an interesting property to be noticed is that the value j™ a:c (0) is independent of the applied magnetic field at 
least as regards the existence of the peak effect in the transport current density (Fig. [7]). Thus, the enhancement of 
the transport current density can be either accelerated or decelerated with the tilt of the applied magnetic field, but 
in general terms, its maximum directly relates to the limitation introduced by the cutting mechanism. Physically, this 
means that the role played by the magnetic anisotropy of the material law may be characterized by the influence of 
the threshold cutting value on the enhancement of the critical current density. This point will be made clearer when 
material laws with different values of x are considered. 
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FIG. 8. (Color online) Magnetic moment (M x , M y ) of the slab as a function of the applied magnetic field component h y (a) for 
the CT1 model. The curves are labelled according to the perpendicular magnetic field component h z o for each. 



3. Magnetic moment of the sample 

Finally, notice that the limitation introduced by the flux cutting mechanism imposes a maximal compression of 
the current density within the sample. Thus, the peak effects both for the transport current density j y (Fig. [7]) and 
for the magnetic moment component M x (Fig. [5J are defined by the instant at which the maximal transport current 
density occurs. Additionally, upon further increasing the longitudinal applied magnetic field component h y (a), the 
profile h x (z) will be forced to decrease from the central sheet {z = 0) towards the external surface (z = a). This 
reversal generates a local distortion of the longitudinal current density j y in a bow tie shape (see the middle row of 
Figs. [5] and Likewise, as soon as the profile j||(0) = j c n) is reached, the magnetic moment M x starts decreasing 
(Fig. [5]) as one can see by comparison of Figs. [7] and [5] 



C. CT^-states: Emergence of the cutting threshold 



The intrinsic interplay between the cutting and depinning mechanisms introduces difficulties when one tries to 
extract the threshold values for a typical experimental arrangement. Nevertheless, as it was argued in the previous 
section, several fingerprints of the actual material law can be identified, i.e.: the longitudinal current density along the 
central and surface superconducting layers, and the behavior of the magnetic moment curves. In order to complement 
this scenario, we have simulated three additional experiments, defined by the so-called CT2, CT3, and CT4 conditions, 
in which the threshold value of the cutting current density is varied relative to the depinning limit (x = 2,3,4 
respectively). 



1. Field and current density penetration profiles 

The behavior of the local electrodynamic quantities is displayed in terms of the aforementioned three stages. In 
order to allow comparison, the main features displayed in Figs. [5] and [5] for the CT1 condition with Ii z q = 0.5 and 
10, are also shown in Figs. 151 and [TU1 under the conditions CT2, CT3 and CT4. Several peculiarities are to be noted. 
To start with, the emergence of negative states for the transport current density close to the external surface of the 
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FIG. 9. (Color online) Profiles of h x [z,h y (a)] and j y [z,h y (a)] for the 1st (top) and 2nd (middle) stages 
of the magnetic dynamics described in CT2, CT3, and CT4 cases (see text) all under a field h z o = 0.5. 
The 3rd stage is only defined for the CT2 case (bottom). h y (a) is as follows. CT2: 1st stage h y (a) = 
0.005,0.050,0.170,0.340,0.500,0.680,0.845,1.0,1.1, 1.3,1.6,1.9,2.2, 2nd stage h y {a) = 3,4,5,6,8,10,20,40, and 3rd stage 
h y (a) = 80, 120, 160, 200, 300, 400, 600, 800, 1000. CT3: 1st stage h y (a) = 0.005, 0.050, 0.170, 0.340, 0.500, 0.680, 0.845, 1.0, 1.5, 
1.8,2.2,2.6,3.0, and 2nd stage h v (a) = 4,5,6,8,15,20,40,70,100,150,200,300,400,600,1000. CT4: 1st 

stage h y (a) = 0.005,0.050,0.170,0.340,0.500,0.680,0.845,1.0,1.5, 1.8,2.2,2.6,3.0,4.0, and 2nd stage h y (a) = 
5, 6, 8, 15, 20, 40, 100, 200, 400, 600, 1000. 



superconducting sample will be more evident either when h z o is reduced and/or the current density anisotropy factor 
X is increased. 



2. Analysis of the current density behavior 

In order to confirm the above interpretation, we show the magnetic dynamics of the longitudinal current density 
j y , the cutting current component j\\ (Fig. ITTj) for the conditions CT2, CT3, and CT4. We have taken a wide set of 
values for the perpendicular field component (h z0 ). On the one hand, as regards the sample's surface, Fig. 1111 shows 
that the longitudinal current density j y (a) does not display significative differences when one has \ > 2 (see also 
Fig. [3] for x = °°)- Thus, the disappearance of the negative current flow along the external superconducting surface 
does not occur despite a very high applied magnetic field has been considered (h y (a) = 1000). On the other hand, it 
is important to notice that the patterns of the parallel current density along the superconducting surface (j\\(a)) are 
almost indistinguishable as soon as the condition x > 2 (CT2) is reached (upper half of Fig. [TTj) . This implies that 
for an accurate picture of the parallel critical current, surface properties do not provide a useful information. 

However, outstandingly, Fig. [11] shows that the threshold value for the cutting current density can be estimated 
from the experimental measurement of the transport current density along the central sheet of the superconducting 
sample. Notice further that, regardless of the experimental conditions (h Z Q, h y {a)) and also for different bandwidths 
X no significant change occurs in the parallel current density around the central sheet of the sample (lower half of 
Fig. HI]). 
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FIG. 10. (Color online) Similar to Fig. [9] but for a perpendicular magnetic field component h z o = 10. The curves 
are labelled as follows. (CT2 - left panel): 1st stage h y (a) = 0.005,0.050,0.170,0.340,0.500,0.680,0.845,1.0,5.0, 
10,15,20,30, and 2nd stage h y (a) = 40,60,80,120,160,200,300,400,600,800,1000. (CT3- middle 

panel): 1st stage h y (a) = 0.005,0.050,0.170,0.340,0.500,0.680,0.845,1.0,5.0, 10,15,20,25,30,40, and 

2nd stage h y (a) = 60,80,120,160,200,300,400,600,800,1000. (CT4 - right panel): 1st stage 

h y {a) = 0.005,0.050,0.170,0.340,0.500,0.680,0.845,1.0,5.0, 10,15,20,35,30,40,60, and 2nd stage h y (a) 
80, 120, 160, 200, 300, 400, 600, 800, 1000. 



3. Magnetic moment of the sample 



The peak effects observed both in the local transport current density j y (0) (lower half of Fig. ITTjl . and on the global 
magnetic moment M x (Fig. I12[) . are predicted to appear subsequent to the maximal longitudinal transport current 
density j™ ax (Eq.pp. 

Furthermore, one additional feature is to be noted: the interval between the instant at which the maximal transport 
current condition is reached (j y (0) = j y nax ), and the instant at which the slope of the magnetic moment M x changes 
sign, could be assumed as transient or stabilization period required for an accurate determination of the value j c m 
when measurements are performed in terms of the applied longitudinal field h y . Apparently, this transient increases 
with the value of the perpendicular field h z o- From this point on, j y (0) may be basically identified with j c ^. 



IV. CONCLUSIONS 



Despite of extensive experimental and theoretical studies about the electrodynamic response of type-II super- 
conductors in longitudinal geometries, much uncertainties remain about the interaction between flux depinning and 
cutting mechanisms, and their influence in such striking observations as the appearance of negative transport current 
flow, the enhancement of the critical transport current density, and the observation of peak effects on the magneti- 
zation curves. In this article, and based on the application of our general critical state theory, we have reproduced 
theoretically the existence of negative flow domains, local and global paramagnetic structures, emergence of peak-like 
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FIG. 11. (Color online) Dynamics of the current density vector for the CT2, CT3 and CT4 models. In the upper panel, we 
show the patterns of transport current along the central (j y (z = 0)) and surface layers (j y (z = a)) of the slab. The curves are 
labelled according to the perpendicular magnetic field component h z o = IE — 7, 0.5, 1, 2, 5, 10 



structures in the longitudinal magnetic moment, as well as the compression of the transport current density for a wide 
number of experimental conditions. 

The longitudinal transport problem in superconducting slab geometry has been studied as follows: we have consid- 
ered a superconducting slab lying at the xy plane and subjected to a transport current density along the y direction. 
The slab is assumed to be penetrated by a uniform vortex array along the z direction, so that the local current 
density along the thickness of the sample is entirely governed by the depinning component Jj_, perpendicular to the 
local magnetic field. Subsequently, a magnetic field source parallel to the transport current direction is switched 
on. Then, the experimental conditions have been modulated through the value of the external magnetic field H v q. 
The dynamical behavior of the transport current density is shown to rely on the interaction between the cutting and 
depinning mechanisms. Moreover, the intensity of the inherent effects has been shown to depend on the perpendicular 
component H z o, being more prominent as this quantity is reduced. By means of our general critical state theory that 
allows to modulate the influence of the different physical effects, we have been able to show that the peak structures 
observed in the magnetization curves and the patterns of the transport current along the central section of a super- 
conducting sample arc both directly associated with the local structure of the vortex lattice. Such dependence may 
become more pronounced as the extrinsic pinning of the material is reduced, in favor of the flux cutting interactions. 
The same conclusion was pointed out from the experimental measurements of Blamire et. al (ReiJ^-) for high critical 
temperature and conventional superconductors. By using a theory that depends on two parameters (x = J c ||/J c j_ 
and n) accounting respectively for the intrinsic material anisotropy and for the smoothness of the Ju(J±) law, here, 
we have quantitatively investigated the influence of the flux cutting mechanism. It has been done by comparing the 
so-called square model (x = 1), the T-state model {\ — > oo), and the double critical state conditions CTx with \ = 2, 
X = 3, and % = 4, all of them with the smoothing index n = 4. 

From our theoretical framework we obtain that the isotropic model (circular region: x = lj n — 1) does neither 
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predict the appearance of negative current patterns nor the peak effects in the magnetic moment curves. However, 
as long as a clear distinction between the depinning and the cutting components of J is allowed, several remarkable 
facts can be explained. 

Going into detail, when the cutting threshold is high (J c ii 3> J c j_ or X ^ 1) the emergence of negative current 
patterns is ensured because unbounded parallel current density allows unconstrained rotations for the flux lines as 
the longitudinal magnetic field increases. Thus, under a range of conditions, the peak effects in the magnetic moment 
and a modulation of the negative surface currents are predicted. 

Concentrating on the local properties within the sample, a clear independence of the field penetration profiles 
relative to the anisotropy level of the material law has been obtained for the partial penetration regime. On the other 
hand, as soon as the full penetration state is reached, noticeable effects of the magnetic anisotropy law are predicted 
both within the central and external layers of the superconducting sample. Thus, according to our results (see Fig |ll[) 
the elusive parallel critical current density parameter J c | can be obtained from the local measurement of the transport 
current density along the central layer of the superconducting sample. 
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